Single-photon excitation of surface plasmon polaritons 
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We provide the quantum mechanical description of the excitation of surface plasmon polaritons 
on metal surfaces by single-photons. An attenuated-reflection setup is described for the quantum 
excitation process in which we find remarkably efficient photon-to-surface plasmon wavepacket- 
transfer. Using a fully quantized treatment of the fields, we introduce the Hamiltonian for their 
interaction and study the quantum statistics during transfer with and without losses in the metal. 
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The emerging field of plasmonics ^ is experiencing 
a considerable increase in interest from researchers in 
many areas of the physical sciences Q ■ Plasmonic-based 
nanophotonic devices in particular have begun to at- 
tract a keen interest from the quantum optics commu- 
nity for their use in quantum information processing 
(QIP) d, [3j S In order to unlock the potential that 
plasmonics at the quantum level can offer, a clear un- 
derstanding of the interplay between single-photons and 
surface plasmon polaritons (SPPs) is of fundamental im- 
portance. Recent studies have focused on various sys- 
tems where SPPs and photons interact 0, IE 0| • How- 
ever, a major obstacle has been the low transfer efficien- 
cies found at the single-photon level [1, Q , thus a com- 
plete quantum description of an efficient transfer process 
is highly desirable. With an extensive understanding of 
photon-SPP coupling in the quantum regime, we can ex- 
pect to open up an array of new applications in QIP, 
based on linear and nonlinear plasmonic effects facilitated 
by strong electromagnetic field confinement [5|, |7| . 

In this work we provide the first quantum description 
of the coupling between single photons and SPPs in a ver- 
satile attenuated-reflection (ATR) setup previously used 
only for classical SPP generation [§] . This is dis- 
tinct from earlier work, such as couplings at rough sur- 
faces [1C| , requiring an entirely different approach. The 
Hamiltonian that we introduce is based on a fully quan- 
tized treatment of both photon and SPP field modes and 
applies to a wide-range of ATR parameters. We find that 
remarkably high quantum efficiencies can be reached for 
photon-to-SPP transfer. We then establish the extent to 
which the excited SPPs preserve the quantum statistics 
of the photons as they travel on realistic metal surfaces. 
Our work provides significant insights into the physics of 
photon-SPP coupling at the quantum level. The methods 
developed are well-suited to other coupling geometries. 

SPPs are highly confined, nonradiative electromag- 
netic excitations associated with electron charge density 
waves propagating along a dielectric-metal interface. In 
Fig- m (a) we show the ATR setup utilized for single- 



photon excitation of SPPs. At various points we will 
introduce the metal as silver only to illustrate our main 
results; the theory developed here fits a far more general 
setting. For SPP excitations, due to the collective na- 
ture of the electron charge density waves, a macroscopic 
picture of the resulting electromagnetic field is appropri- 
ate [lol |. Upon quantization, SPPs therefore correspond 
to bosonic modes. The quantized vector potential in the 
continuum limit for SPPs propagating along an air-metal 
interface at z = in the x direction^ as shown on the 
right of Fig. [T] (a), is given by [13, IIH kspp{r,t) cx 
duj{M{io)L)-^/'^[(j){r,uj)e-^'^%io) + h.c]. The disper- 
sion relation is uj"^ — c^lP + / tm with e,„ the permit- 
tivity of the metal, M{uj) is a frequency dependent nor- 
malization [10] and L is the profile-width flU. The 5(a;)'s 
(6^(ci;)'s) correspond to annihilation (creation) operators 
which obey commutation relations [6(w), w {oj')] — 6{uj — 
uj'). The mode functions are given by 4>{r,u!) = [{ix — 
ki/u)e-'''-&{z) + {i± + fcz/^^o)e'''''=^?(-^)]e*■^ where the 
wavevector k — /cx, 'd{z) is the Heaviside step function 
and the decay of the SPP into the metal (air) is parame- 
terized by v'^ — fc^ — EmW^/c^ (z/q = k'^ — uj"^ / c^) . For pho- 
tons propagating in air in the k' direction (k' = sin 0x -|- 
cos6'z, as shown on the left of Fig. [T] (a)), we have 11| 
Ap(r,i) cx /o°°dw(wA)-i/2[e''='(k'-r)e-«-ta(^) + f^.c]. 
The dispersion relation is w = cfc', A is the beam cross- 
section and [a{uj) , a"^ {lu')] = S{uj — lu'). Here, the SPPs 
and photons are transverse magnetic modes. At the 
single-photon level only small intensities of the photon 
field are involved and any nonlinear terms in the photon- 
SPP couphng can be sufficiently neglected [l3|. We are 
thus led to the following natural linear coupling Hamil- 
tonian for the entire system shown in Fig. [T] (a). 



Us = 



dujhuja) {uj)a{uj) -\- i dojfujjW {uj)b{ijj) (1) 



duj[g{uj)a) {ijj)b{uj) — g* {ijj)b\uj)a{uj)]. 



The first and second terms are the photon and SPP fields' 
free-energy respectively. The last term, which we denote 
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FIG. 1: Single-photon excitation of SPPs using attenuated- 
reflection, (a): A photon wavepacket is injected into the 
system at a specific angle 9, with a prism mediating an inter- 
action between the photon and SPP modes. The minimum 
prism size is diffraction-limited, (b): Dispersion relations for 
the photon (shaded region) and SPP (curve). The prism en- 
ables mode-matching, (c): Two ATR excitation geometries; 
(i) Otto and (ii) Kretschmann-Raether, see text for details, 
(d) : Transfer process for the photon and SPP mode operators. 

as Hint, describes interactions between the two fields, 
where the couphng g(w) is a function of the system pa- 
rameters for a given ATR geometry. In Fig. [T] (b) we 
show the dispersion relations for SPPs and photons. As 
an example, we choose em — l^'^p/'^^ + ^^mi where ujp — 
1.402 X lO^^rad/s is the plasma frequency for silver and 
de^^ = 29u}^ /ujp is a background correction term Ne- 
glecting Ohmic losses, over the uj range of the SPP modes, 
k produces a curve which approaches the surface plasma 
frequency u>sp, where e™ = —1. On the other hand, the 
X component of k' for photons in air incident at angle 9 
covers the shaded region of Fig. [1] (b) . Two ATR geome- 
tries that can provide the necessary mode-matching for 
coupling photons to SPPs are shown in Fig. [T] (c), de- 
noted as (i) Otto (O) and (ii) Kretschmann-Raether 
(KR) [9]. Both consist of a prism in layer I with per- 
mittivity ei. The O (KR) geometry has air in layers II 
and IV (III and IV) with £2 = £4 = 1 (£3 = £4 = 1) and 
metal in layer III (II) with £3 = £„ (£2 = £m). In both, 
SPPs are excited on the II/III interface, with z ^ z — d 
in (j){r,oj). For the KR geometry: 1/ ^ vq. The thick- 
ness I is assumed to be far larger than the decay of the 
SPP into the metal (air), i.e. I ^ ('^o'^)' making the 
effects of layer IV negligible. 

With the ATR setup introduced, we can now formulate 
the photon-SPP coupling model of Eq. ([1]) . For each lo 
in both geometries the couphng can be described by a 
transfer matrix T{lj) in the Heisenberg picture as 

The transfer process is depicted in Fig. [T] (d), where the 
commutation relations of the quantum operators a{Lo) 
and &(w) define the structure of T(a;), while its coeffi- 



cients (|a(a;)p-|-|/3(w)p = l,Va;) are determined from the 
overlap of system mo defunct ions. By solving Maxwell's 
equations across the first three layers shown in Fig. [T] (c) , 
one finds the modefunctions of the field in layers II 
and III: il^ir.ij) = {[{ipie-^'^'' + i^2e''"^)x + (ipze'^'^'' + 
V?4e'^""^)z;]i?(z)i?(d ~ z) + [(^se^'^^^x + <^^e-'^^'-'L\d{z - 
^^|gjKx_ jjg]-g^ i^jjg (p.'g g^]-g constants related by bound- 
ary conditions at the interfaces, 7,; = [k^ — tiUj'^ / c^Y^"^ 
and the dispersion relation k = y/ei{u)/c)sa\6. The O 
(KR) geometry has £2 = 1 (£m) and £3 = £m (1) with 
Em = 1— Wp/(<^^+«wr)-|-(5£m, which now includes a damp- 
ing factor r for the metal and a complex correction term 
The complete modefunctions for the three-layer 
(3L) system are: '^{t.lo) — r'i(j{r,uj)-d{—z) + Ttl;{r,uj)'d{z) 
where r and t (|rp + jrp = 1) are obtained from Fres- 
nel's relations at the boundaries. However, the ip{r,uj) 
are not involved in the coupling due to mode-matching; 
they always have a real component of their wavevector 
in z. On the other hand, mode-matching can be satisfied 
between the two-layer (2L) modefunctions ip{r, uj) and 
the SPP modefunctions by fixing the angle 9 correctly. 
For instance, by setting k = k the dispersion lines cross 
at 6* = sin~^[£m/(ei(l + Em))]^^^ in both geometries. In 
Fig. [H (b) we show this for a particular angle 9 = 85° 
(2L line). The inset shows that mode-matching over the 
entire range of uj can be achieved, for example, using a 
prism with £1 — 1.51 and silver with T = 6.25 x lO-'^^rad/s 
and Ssm = Sem+i^^lm where Se^ — 0.22 fl^. This range 
is important for excitation with a photon wavepacket of 
finite width, as we show later, and is not possible in other 
excitation schemes such as the grating-type coupler. 

In Fig. [1] (d) the bin/outi'^) operators are associated 
with the in/out SPP modefunctions 0(r,cj) (bin{uj) = 
b{uj)) and ciin/outi^) with the in/out 3L modefunc- 
tions ^'(r, w). For negligible loss on entry into the 
prism medium, we can assume the operator relation 
ain{uj) = a{uj). We then have (3*{uj) = —T[S{LJ — Lj')6{k — 
^)Jdz{J^~'^\uj)^{r,u)r ■ {J^^'/\uj')i^{r,uj'))] ^ ^. 
Several factors permit the use of the mode overlap in the 
value of f3*{Lj). First, we assume that the SPP modes 
experience negligible damping during the excitation pro- 
cess, Im(£„i) w 0, imposing damping effects subsequently 
as the SPP propagates. Second, the SPP is assumed to 
exit the prism region on a time-scale such that mode- 
matching conditions are broken almost immediately af- 
ter excitation. This can be achieved by adjusting the 
excitation point Third, as the SPP modefunctions 
exist in the region z £ (— oo,cx)), d must be chosen 
such that their decay into the prism is negligible, allow- 
ing it to be neglected from their definitions. To check 
an acceptable range of d we define a penetration-factor 
p = 2/uod {2/ud) for the O (KR) geometry and consider 
the SPP modes as good approximations for V <1, where 
|(/)(r, uj)\'^ at z = is less than 2% its maximum value. In 
Figs. (a) and (b) we use the example of silver to show 



3 



V over a range of u> and d for the two ATR geometries. 

With the above considerations, we can now determine 
the couphng g{uj). In order to connect (}*{u)) and 
we set a = cos 6 and (3 = e**sine (9 e [0,7r/2], 
$ S [0, 27r]) for each uj and parameter set {d, £^,6'} 
for a given geometry. The corresponding Hamiltonian 
in the Schrodinger picture is Tiint from Eq. ([1]) with 
g{io) = e'*('^)e(w) Q. In Figs. [2] (c) and (d) we use 
silver to plot a rescaled couphng, \g{oj)\ — ^\g{uj)\, for 
the two ATR geometries, such that = 1 (0) corre- 

sponds to a unit (zero) transfer probability of a photon 
to an SPP. Figs. [2] (a) and (b) show the optimal 
for both geometries satisfy V ^ \. In Figs. [2] (e) and 
(f) we plot these optimal values and the value of d at 
which they occur. The optimal |5(w)| in both geome- 
tries rises for increasing w, reaching an apex, then drops 
sharply as uj tends toward oj^p- This behavior is due 
to a dominance of the value for r in P*{lo) at large to, 
which decreases rapidly due to boundary conditions and 
the large required for mode-matching. Such excellent 
coupling values have been found classically [%, 9], how- 
ever, this is the first time a rigorous quantum mechanical 
treatment has been achieved, making it possible for us 
to determine correctly the quantum efficiency of single- 
photon excitation. The coupling cannot be deduced 
from a classical model of the system. 

So far we have focused on single modes of the system. 
However, it is important to consider the transfer of a 
photon wavepacket, such as in an experiment, to an SPP 
wavepacket state. An n-photon wavepacket state is given 
by |nj) = (n!)-i/2(at)n|o)^ ^here d\ = ^ duj£,{uj)a) {lo) 

with J da;|^(cL')p = 1 [lH. For simplicity we take a 
Gaussian profile for a wavepacket produced at time 
to = with bandwidth Atj — 2a \/ 2 log 2 and central fre- 
quency ojQ. We allow to e (—00,00) as Aw <C oj. For 
wavepacket-transfer with negligible deformation^ each lo 
must have approximately the same g{uj) and 9. This 
is satisfied given a small enough bandwidth with slowly 
varying g{uj) and 9. For example, using silver with 
AA = lOnm, one finds that A0 rises exponentially from 
0.004° at 1 X lO^^rad/s to 14.61° at 5 x lO^^rad/s, which 
can be attributed to the dependence of on w (see inset of 
Fig. [11(b)). The couplings show a similar behavior, with 
|A(7(w)| = 0.01 for d optimizing |5(wo)| at 1 x lO^^rad/s, 
rising to \l^g{uj)\ = 0.2 (0.04) at 5 x W^va<l/s for the O 
(KR) geometry. The large difference is due to the sharper 
drop in |f/(w)| for O at high lo. Significant wavepacket 
deformation can be avoided by operating at low uj, al- 
though at the expense of the coupling. In general, a 
narrow bandwidth will provide access to larger couplings 
with negligible deformation. 

Finally we turn our attention to damping in the metal 
as the excited SPP propagates. For the couphng, the 
approximation Im(e„i) « was made for the excited 
SPP. However as it travels, finite conductivity of the 
metal and surface roughness result in heating and radia- 
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FIG. 2; Photon-SPP coupling; The left (right) column corre- 
sponds to the O (KR) geometry. Panels (a) and (b) depict 
the penetration-factor V. Panels (c) and (d) show the behav- 
ior of the coupling \g{ijj)\ = §|(;([i;)|. The optimal values of 
|g(tj)| (dashed line) are displayed in panels (e) and (f) along 
with the values of d (solid line) at which they occur. 



five losses respectively [ij; for a reasonably smooth sur- 
face, thermal loss is the main source of damping. While 
a quantization of the decayed SPP modes can be per- 
formed, a mathematically equivalent and simpler model 
is the method of arrays shown in Fig. [3] (a) [11]. Here 
we introduce a bath of field modes, described by op- 
erators Ci(w) {i = l,..,iV) separated by Ax, interact- 
ing with the SPP wavepacket as it propagates. In the 
limit N ^ 00 and Ax 0, the SPP operator becomes 
bZtito) = e^^4„„,(c^) + iV2^/;dx'e*^(--')c(c^,x'), 
with Ci{Lij) — > V Axc{u), x') and 5ij — s- Ax5{x' — x"). 
The array coefficients are chosen such that the bath 
modes induce a change in the SPP wavevector k match- 
ing that of the complex i.e. k ^ K = (w/c)[em/(l + 
^m)Y^'^ = k + iK{uj), where 2k{lj) is the loss per 
unit length of propagation. We then set the relations 
{c{lu,x')) = (ct(a;,x')) = {c^ iLU,x')c{LJ,x")) = for 
the bath modes at room temperature and the frequen- 
cies considered [lH. We assume that the excited SPP 
wavepacket with wq has a narrow enough bandwidth 
such that k(lu) sa k^luq) — kq and k w k^LUo) + (lu — 

u;o)vq'^{ujo), with Vq^{ljo) = ^l^j^L^t^o- The flux of 
SPPs at point x along the metal surface is then sim- 
ply foutit) = {b^^t)b^ut(t)) = e-^-o'^{bl,{tn)bout{tR)), 
where ta = t — xvq^{luo). For an initial SPP wavepacket 

with n excitations, (SLtiip)L«t(^fl)) = '^ll^(*fl)P- A de- 
tector with efficiency fj, [111 . Il7l | operating for time period 

[xv'ij^{ljo) — 1/cr, xv'^^Iujq) + 1/a] would measure a mean 
SPP-count of (to) = J dt foutit) = //ne-2«o^. The 
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FIG. 3: Damping model for propagating SPPs. (a): Bath 
modes interact with the SPP mode as it travels along the 
metal. Panels (b) and (c) show the normalized expected 
mean SPP-count (me)/n at point x for an n photon-to-SPP 
wavepacket transfer via the O and KR geometries respec- 
tively. 



We have provided the first quantum description of the 
Photon-to-SPP transfer process for ATR excitation. Re- 
markably good quantum couplings over a wide-range of 
frequencies were found. We also examined the extent 
to which the excited SPPs preserve quantum statistical 
properties. The techniques developed here provide key 
insights into the formulation of quantum descriptions for 
the photonic excitation of SPPs. This work can therefore 
be seen as an important starting point for future research 
into the design of new quantum plasmonic devices for ap- 
plications based at the nanoscale, such as SPP-enhanced 
nonlinear photon interactions and SPP-assisted photonic 
quantum networking and processing. 
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Photon-to-SPP transfer process must however be incor- 
porated for determining the expected mean SPP-count 
(me) from an incident n photon wavepacket. The en- 
tire process is analogous to an inefficient detection prob- 
lem [13] and we have (me) — ^|/3(a;o)pne~^''°^. In 
Figs. [3] (b) and (c) we show (m^jn for the ATR ge- 
ometries, where [i, = 0.65 is chosen as an example of 
non-ideal signal extraction using, 6.17. a prism and pho- 
todetector. The detection of many identical excitations 
from a set rate of single photons would be required to 
determine (me). 

While the quantum observable (me) matches well the 
behavior of its classical counterpart, the field intensity 
/ it is not sufficient in an experiment to show that 
the SPPs are quantum excitations. We now consider an- 
other observable, the zero time-delay second-order quan- 
tum coherence function (;(^)(0) [ll| at a fixed position, 
defined as g(2)(o) = {-.Pit):) j {: i{t):f . Here, / is 
the intensity of the quantized field operator, : : de- 
notes normal-ordering and the expectation value is taken 
over the initial state of the field. For a classical field 
1 < ff^^''(0) < 00. On the other hand, for an inci- 
dent n photon wavepacket g*-^^(0) — {m{m — l))/(m)^, 
where (m) = n J^*"'"^ '^+'it''+''>i2 

rt+T 



n 



mt^dt'm')??, giving g^^^m = 1 - l/n. 



and {m{m — 1)) = n{n — 
(2)rn^ = 1-1 /n. This al- 
ways lies in the classically forbidden region (7^^'(0) < 1. 
The value of ^^^•'(O) for an excited SPP wavepacket at 
point X can be found by recognizing that the photon- 
to-SPP transfer and SPP propagation stages constitute 
an array of lossy beamsplitters 17|. At a beamsplit- 
ter with loss coefficient rj^^"^, the quantum observables 
(m) — !■ 77(m) and {m{m — 1)) — > if'{m{ra — 1)). Thus, 
the individual losses accumulated cancel, leaving ^^^^(O) 
surprisingly unaffected. A Hanbury-Brown Twiss type 
experiment [18|] could be used to measure g^^^(O). 
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